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ABSTRACT 

We construct a holographic map between asymptotically AdS^ x solutions of lOd 
supergravity and vacuum expectation values of gauge invariant operators of the dual QFT. 
The ingredients that enter in the construction are (i) gauge invariant variables so that the 
KK reduction is independent of any choice of gauge fixing; (ii) the non-linear KK reduction 
map from 10 to 5 dimensions (constructed perturbatively in the number of fields); (iii) 
application of holographic renormalization. A non-trivial role in the last step is played by 
extremal couplings. This map allows one to reliably compute vevs of operators dual to 
any KK fields. As an application we consider a Coulomb branch solution and compute 
the first two non-trivial vevs, involving operators of dimension 2 and 4, and reproduce the 
field theory result, in agreement with non-renormalization theorems. This constitutes the 
first quantitative test of the gravity/gauge theory duality away from the conformal point 
involving a vev of an operator dual to a KK field (which is not one of the gauged supergravity 
fields). 
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1 Introduction 

Gravity/gauge theory dualities relate string theory on spacetimes that asymptote to AdSm x 
X, where X is a compact manifold and gauge theory residing on the conformal boundary 
of the AdS part of the geometry. In the initial work [1] the dual theory was a conformal 
field theory (CFT) and the bulk spacetime AdSm x X (rather than asymptotic to it), but 
it was soon recognized that the duality can be extended to describe quantum field theories 
that can be obtained from the CFT by either adding new terms in the action or considering 
vacua where the conformal symmetry is spontaneously broken. Both of these cases are 
described gravitationally by a solution that is asymptotic to AdSm x X. 

Despite much work however basic questions still remain. One such question that will be 
the subject of this paper is: 

Given a ten dimensional solution that is asymptotic to AdSm x X how does one compute 
the vacuum expectation values of gauge invariant operators ? 

Roughly speaking vevs of chiral primary operators should appear in the radial expansion 
of the bulk solution. However making this precise proves to be a lot more subtle that one 
might have anticipated, and even qualitative features are not reproduced correctly by naive 
methods. The answer to this question should follow from the basic AdS/CFT dictionary 
[2, 3]. This is indeed the case but in order to implement the idea one has to sharpen existing 
methods and overcome several technical issues. 

To illustrate the issues involved it is instructive to consider a simple example where 
the physics of the solution is well understood. A class of such examples is provided by 
multicenter D3-brane solutions, which in the near-horizon limit correspond to the Coulomb 
branch of A/" = 4 SYM [4] . These examples are particularly interesting because the QFT 
vevs are protected by a non-renormalization theorem, and the gravitational results must 
therefore agree exactly with those computed at weak coupling. The metric is of the well- 
known form 



where H is a harmonic function in transverse directions. For a distribution a{y) of D3 
branes, the harmonic function reads 



where in the last equality we expanded in = |a;_Lp, are spherical harmonics and Qj^ 
are numerical constants that depend on the distribution a{y). Inserting this in (1.1) and 



c/s^ = H{x^)-^/'^dxl + H{x^Y''^dx\ 



(1.1) 




(1.2) 
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expanding in r results in a metric that is asymptotically AdS^ x S^. The QFT vevs should be 
encoded in the asymptotics and the purpose of this work is to show how to unambiguously 
extract this information. 

The solutions under discussion are special in that they are uniquely determined in terms 
of a harmonic function. Furthermore, the spherical harmonics appearing in (1.2) are in 1-1 
correspondence with the chiral primary operators of = 4 SYM and the radial power at 
which they appear is the correct power for their coefficients to correspond to the vevs of 
the dual operator. This led [5] to propose that these coefficients are proportional to vevs 
of the corresponding operators. Although this is a well motivated proposal, it is not clear 
how one would generalize it to the general case where the solution is not determined by a 
harmonic function. Even for the case at hand there are several open questions. For instance, 
inserting the harmonic function in the metric leads to terms involving powers of spherical 
harmonics whose meaning is not clear and in general there is also a dependence on the 
radial coordinate used to perform the asymptotic expansion. The simplifications special to 
such cases (i.e. when the solution is determined by a harmonic function) will be discussed 
in a separate publication [6]. In this paper we strive for generality, so our starting point 
will be general asymptotically AdS^ x metrics and we will only use the CB solutions in 
order to illustrate the general procedure. 

Recall that the basic dictionary of the gravity /gauge theory duality [2, 3] states that 
(i) there is a bulk field corresponding to each gauge invariant operator and (ii) the string 
partition function with bulk fields satisfying appropriate boundary conditions is equal to 
the generating functional of QFT correlators with the boundary conditions playing the 
role of sources. In particular, one could compute vevs (in the large N and large 't Hooft 
coupling limit) by differentiating the supergravity on-shell action once w.r.t. sources. In 
implementing this procedure however one finds several obstacles. 

First, the relation in (ii) should be understood as a "bare relation" as both sides di- 
vergence. To make the procedure well-defined one must renormalize. This is a standard 
procedure on the field theory side. On the gravitational side, the corresponding procedure, 
denoted holographic renormalization, was developed in a series of papers [7, 8, 9, 10] (see 
also [11]-[16] for related work and [17] for a review)^. After renormalization is done, the 
one point functions can be computed in all generality. The answer relates the one point 
function to certain coefficients in the asymptotic expansion of the bulk fields. So given any 

^The starting point in the analysis in these papers was the lower dimensional AdS gravity obtained by 
reducing the original theory over the compact space X. A discussion that starts from higher dimensions can 
be found in [18, 19]. 
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solution of 5d gravity coupled to matter one could read off the vevs of the dual operators 
by looking at the asymptotics. 

We would like to emphasize that the procedure of renormalization is essential for cor- 
rectly extracting the vevs. To give an example where a naive prescription fails, consider 
the case of the CB solution corresponding to a distribution of D3 branes on a disc of radius 
I. The 5d metric obtained by reducing the lOd solution over the sphere has the following 
asymptotics, 



A naive prescription for reading off vevs that is often quoted in the literature is that the vev 
of an operator can be obtained, up to a (non-zero) numerical constant, from the normalizable 
mode of the corresponding bulk field. The bulk metric is dual to the stress energy tensor 
so one would be tempted to identify the coefficient of the term with the vev of the stress 
energy tensor. This is clearly incorrect since that would imply non-zero vacuum energy 
for the dual theory but the solution is supcrsyninictric so the vacuum energy should be 
equal to zero. Indeed the 1-point function extracted using holographic renormalization 
[8, 9] contains additional terms (sec (5.11) below) and taking those into account one finds 
the expected result, zero. Such subtleties are present in all cases, including that of scalar 
fields. For example, for deformation fiows the vevs of all operators should be equal to zero, 
but there are examples where the above naive prescription leads to non-zero values. Again 
the correct 1-point functions include additional terms so that the total result is zero (see 



An analysis that starts from the lower dimensional gauged supergravity is suflficient if one 
is only interested in computing vevs for operators dual to fields of the gauged supergravity. 
There is however an infinity of other (half supersymmetric) gauge invariant operators which 
would then ab initio be excluded from the analysis. These operators are dual to massive 
KK fields. The map between KK fields and gauge invariant operators was worked out in [3] 
(and subsequent papers) using the computation of the KK spectrum of AdS^ x in [21]. 
In this paper a linearized analysis around AdS^ x was performed. This analysis provides 
an explicit map (in a specific gauge) between linearized solutions of the ten dimensional 
equations of motion and linearized solutions of the dimensionally reduced five dimensional 
equations. To compute the vevs however we need to know the map at the non-linear level. 
To illustrate this, let s'' be ten-dimensional fields and let S'' be the corresponding five 




(1.3) 



[8, 9, 20]). 
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dimensional fields. In general the reduction map will be non-linear and takes the form 

= + E {Jklms's"' + LkimD^s'D^^s^ + 0[sf) (1.4) 

Im 

where Jtim and Lkim are numerical coefficients and we retain only terms quadratic in the 
fields. (We also suppress contributions on the right hand side from other scalar and non- 
scalar fields since they are not necessary to illustrate our point.) If is dual to an operator 
of dimension k then we would need to extract the coefficient of order to determine the 
operator's vev. Clearly, quadratic terms with I + m = k will also contribute at the same 
order and therefore STich non-linear terms in the KK map cannot be ignored. Similarly 
cubic and higher order in fluctuation terms that are of order z*^ will also contribute, along 
with non-linear contributions involving other supergravity fields. 

So to read off the vevs we need to understand the KK reduction map at the non- 
linear level. However, if we are interested in the vev of an operator of a given dimension, 
only certain non-linear terms need to be computed, namely the ones that could possibly 
contribute to the vev. For instance, if we are interested in computing the vev of an operator 
of dimension 4, we would only need to keep terms quadratic in the fields dual to operators 
of dimension 2. This in effect truncates the reduction to a finite number of fields. This 
should be contrasted with the issue of consistent truncation. When the latter is possible 
one keeps only the "massless" KK modes in the reduction. In our case we keep massive KK 
fields as well. However, only a finite subset of them contribute to the asymptotics up to a 
given order. 

Another issue is that of gauge fixing. The analysis in [21] was done in the de Bon- 
der gauge. Generically however a given supergravity solution will not be in this gauge 
and finding the coordinate transformation that would bring the solution to this gauge is a 
complicated task. To deal with this issue we will instead develop a "gauge invariant KK 
reduction". Instead of fixing the gauge, we combine the fields in gauge invariant combi- 
nations. This can be done systematically in an expansion in the number of fields. Having 
worked out these combinations, one can relax the de Bonder gauge condition by simply 
replacing every field by its gauge invariant generalization in all results obtained in a specific 
gauge. 

To summarize, we argue that in order to compute the vevs we need to obtain the non- 
linear KK map in terms of gauge invariant variables to appropriate order in the number of 
fields. This procedure results in five dimensional field equations and an explicit map between 
lOd solutions and solutions of these 5d equations. The 5d equations can be integrated into 
a 5d action and from here one can obtain the 1-point functions following the procedure of 
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holographic renormahzation. There is however an additional subtlety. In some cases the five 
dimensional equations contain no couplings between certain fields but boundary interactions 
exist [22]. These boundary couplings are in fact responsible for extremal n-point functions, 
namely correlators involving operators of dimensions {ki=k2+ ■ ■ ■ +kn,k2, ■ ■ ■ ,kn}. One 
must take into account these additional boundary terms when working out the holographic 
1-point functions. 

Combining the non-linear, gauge invariant KK reduction map with the holographic 1- 
point functions we finally arrive at a well defined map between the asymptotics of a lOd 
solution and vevs of gauge invariant operators. 

This paper is organized as follows. In the next section we discuss the Coulomb branch 
of A/" = 4 SYM. We focus on a specific case where the vevs are uniformly distributed on a 
disc and compute all vevs of gauge invariant operators. The challenge for the gravity/gauge 
theory duality is to reproduce exactly these vevs holographically. In sections 3, 4 and 5 we 
build the holographic map. In section 3 we construct gauge invariant variables; in section 
4 we work out the KK map to second order in the fields and in section 5 we derive the 
holographic 1-point functions. In section 6 we discuss the supergravity solution dual to 
the CB state discussed in section 2 and use the map developed in sections 3, 4, 5 in order 
to compute the first two non-trivial vevs and find perfect agreement with field theory! 
We conclude in section 7 with a discussion of our results. Several technical details are 
relegated to appendices A, B and C. In appendix A we discuss the harmonic expansion of the 
antisymmetric gauge field; in appendix B we summarize and develop several results about 
spherical harmonics with SO (4) symmetry and in appendix C we discuss the computation 
of the field equations to second order in fluctuations. 

2 Af = 4 SYM on the Coulomb branch 

J\f = 4 SYM contains 6 scalar fields X^'^ in the adjoint representation of the gauge group 
that we take to be SU{N). The Coulomb branch of = 4 SYM corresponds to giving a 
vacuum expectation value (vev) to the scalars subject to the condition X'^^] = 0. Upon 
diagonalizing the scalar fields the moduli space is parametrized by the 6(iV — 1) eigenvalues 
of vevs (modulo the Weyl group). In the large N limit we can approximate the eigenvalues 
by a continuous distribution. Notice that the Coulomb branch still preserves N = 4 su- 
persymmetry but the conformal symmetry and the 16 superconformal supersymmetries are 
broken. This implies that the vevs are protected from acquiring quantum corrections, as we 
explain at the end of this section. So this example represents an ideal case for a precision 
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test of the gravity/gauge theory correspondence in a non-conformal setting. It has been 
known since [4] that there is a one to one correspondence between the Coulomb branch 
of = 4 SYM and multicenter D3 brane solutions. Since the vevs do not renormalize, 
however, one should be able to establish a strong result, namely the exact values of vevs 
should be reproducible by a gravitational computation. 

We will consider in this paper the specific case of a Tiniform distribution of eigenvalues 
of and X'^ on a disc of radius o and vanishing vev for the remaining scalars, {X^) = 
{X^} = = = 0. Let 

X^=pcos0, X^ = /9sin0 (2.1) 

To leading order in the large A'^ limit we may represent the eigenvalues by a uniform con- 
tinuous distribution, 

a(p,<^) = ^. (2.2) 
Tra 

Notice that the configuration corresponding to this continuous distribution preserves an 
SO {A) X SO {2) symmetry of SO (6). 

To compare with supergravity we would like to parametrize the moduli space by vevs 
of composite operators. We consider the following chiral primaries (CPOs) of A/" = 4 SYM 

0'^=J\fi,C',l,.TriX^^...X^^), (2.3) 

where Mi^ is a normalization factor and C^^ is a totally symmetric traceless rank k tensor 
of SO{6) which is normalized such that {C^^C^^} = Cj....^Cj....^ = S^^^\ The SYM action 
is normalized such that the relevant propagators are 

^li-my)) = Tirfe^' (2-4) 

' (27r)^ \x — y\ 
where a, b are color indices. 

The cases of interest are the operators which are singlets under the decomposition of 
iS'0(6) into 50(2) x 5*0(4) since non-singlet operators have zero vev. These operators can 

be obtained from the explicit expression of scalar harmonics in appendix B.l by suitably 
replacing x*^ by X^^ , compare (2.3) and (B.ll). The result for the singlets is 

(-)" r^JY^^ ^nJ n\(n + m + l 

I m=0 



where B? = Ei=i(^*)^- The explicit expressions for the lowest dimension operators are 
thus 

= A/2^1V(3p2-i?2). (2.6) 



©4 = M4^Tt(Wp* -Sp'^R^ + R'^y, (2.7) 
4v5 

= Ar6^Tr(35p^-45p^i?2 + 15/92i?^-i?6). (2.8) 
To compute the vev of these operators we now use 

(TV(/)) = £dppj\MP,<t^)p' = J^f"^ (2-9) 

and the identity 



H" E(-r " =(»+i). (2.10) 



im=0 




to arrive at 



2n 



(02n\ ^ ^' ^n^ (2.11) 

This result was derived by a tree-level computation. However it remains uncorrected 
both perturbatively and non-perturbatively. A quantum correction to the vev of the scalars 
X' would result from a non-vanishing tadpole contribution and this would induce a cor- 
rection to the effective potential. However, there are no perturbative or non-perburtative 
quantum corrections to the low energy (2-derivative) effective action of iV = 4 SYM [23, 24] 
so the vevs of the scalars are not corrected. The only remaining issue is operator mixing. 
Indeed, chiral primary operators mix with certain multi-trace operators. However, this is a 
subleading effect^ in 1/N and we are considering the leading behavior. It follows that the 
operators have the same vev (2.11) at strong coupling. The challenge for the AdS/CFT 
correspondence is to reproduce these vevs. 



3 KK reduction with gauge invariant variables 

The IIB SUGRA field equations^ for the metric and 5-form field strength are given by: 

Rmn = \fmpqrsFm^^^^ , F = *F. (3.1) 

These equations admit an AdS^ x solution 

dz"^ 1 

dsl = -^ + -^dx'^^+de^ + sm'^ednl + cos'^ed(p'^ (3.2) 

FfivpcTT ~ ^ixvpuTi Fgjy^^ = ^abcde 

^Thc only exception is the case of extremal operators where the mixing with multitrace operators is not 

subleading [22] . In this paper it was argued that the supergravity fields are dual to the single trace operators 

so these are the relevant operators to consider. 

^The field strength differs by a factor of 4 from the conventions in [25]. Index conventions: M,N, ... are 

lOd indices, /i, f, ... are AdSs indices, a,b, ... are indices, x denotes AdS coordinates and y coordinates. 
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We will consider here solutions that are deformations of AdS^ x such that 



9MN = 9mn + hMN, 
FmNPQR = FmNPQR + fMNPQR- 

The fluctuations can be expanded in harmonics: 



(3.3) 



and 



h{ab)ix,y) 



f^uparix, y) 
fa/xupaix, y) 
fabixi/p {.X,y^ 
fabcfiu {X-iy) 

fabcdfj,{x J 2/) 
fabcde{x,y) 



Y,h%[x)Y'^{y) 

E('^S) (^)^A) (y) + ^ll) i^)D(aY;f (y) + </.J) {x)D^aDb)Y'^ (y)) 
J2^'^{x)Y'Hy) (3.4) 



Y^{bXAx)DaY'^{y) + 4Z)[^6^;^j(x)y/«(y)) 
^(3D[,6^-(x)F[;-(y) - 2b^;^^{x)D^^Y;fiy)) 
Y,mA](^)^ab/'DaY,'^{y) + 36^i°(x)D[„y/4°(y)) 



Y,iDf.bll^{x)eabc/DeY'^{y) + (A^^ - 4)6^^(x)e,fe,/y/Hy)) 
Y,bll){x)A''eabcdeY''{y) 



(3.5) 



Numerical constants in these expressions are inserted so as to match with the conven- 
tions of [21], see appendix A. Parentheses denote a symmetric traceless combination (i.e. 

A(^ab) = l/2{Aab+Aba)-l/5gabAl). F^i,y/5^y(^i^'* andYj^i^ denote scalar, vector and tensor 
harmonics whilst A^^ and A^"' arc the cigcnvahics of the scalar and vector harmonics under 
(minus) the d'Alembertian. The subscripts t, v and s denote whether the field is associated 
with tensor, vector or scalar harmonics respectively, whilst the superscript of the harmonic 
label In derives from the number of components n of the harmonic. 

Not all fluctuations are independent however. Some of the modes are diffeomorphic to 
each other or to the background solution, i.e. certain SHmn and SfMNPQR are generated 
by a coordinate transformation, 

(3.6) 



These, up to terms linear in fluctuations, are given by 

ShMN = {Dm^n + Dn^m) + {Du^^hpN + Dn^^Hmp + DphuN); (3.7) 

S fMNPQR = ^D[m^^FnPQR]S + (5-D[MC'^/ArPQil]5 + ^^Ds fMNPQR)- 



The gauge parameter ^ (x, y) can be expanded in harmonics as 

U^,y) = J2^jl{x)Y''{y); (3.8) 

In much of the previous Htcraturc this issue was dealt with by imposing a gauge fixing 
condition, most notably the de Donder-Lorentz gauge fixing condition 

L>"/t(„b) = D^ha^ = 0. (3.9) 

This amounts to setting to zero the coefficients 5^^^^,^^^^,^^^^ (as can be easily seen by 
inserting (3.4) in (3.9)). Although this gauge is a very convenient choice for deriving the 
spectrum, it is not very well suited for holography since generically solutions will not be 
in that gauge. For this reason instead of gauge fixing this symmetry we will derive gauge 
invariant combinations of fluctuations. This will allow us to switch easily between different 
gauges. 

3.1 Gauge invciriance at linear order 

Wc first discuss gauge invariance at leading order, i.e. wc consider the fluctuation inde- 
pendent terms in (3.7). These transformations map the fluctuations to the background 
solution. Under these transformations the coefficients in (3.4) transform as 

S4>'^)=0, H%)='ii[ly ^^i)=Kv '^^'^ = 2A'^^fe- (3-10) 

It follows that is gauge invariant to this order and for the rest of fields we can construct 
the following gauge invariant combinations 



- 7^'--A'^4>\l, (3.11) 
1 

h\L = h{l - DaBjK - D^Bj\ , Ii ^ 0. 



Bj\ = Bj\ --Da(f>i' 



where we define 



Note that the last formula in (3.11) is only valid for Ii^O (the fields B^^l^^ and (f)^^^-^ exists 
only for Ii > 0, since F° is a constant). For Ii = 0, h^^ is a deformation of the background 
metric and from (3.10) we see that it indeed transforms as a metric. 
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Similarly, the leading term in the 5-form transformation implies for the coefficients in 
the harmonic expansion the following transformations 

Sh^i — f -fi Sh^i — f Tfli Sh^5 — \ n 

so that the gauge invariant combinations are 

= ^5)-^^^) (3-14) 
ft^i — ft^i — f 

bj\ = bi\ -^TT- (3.15) 

{v)ij, {v)p, 2{A.^^ - 4) (^) 

and the fields h^^^^p and h^^° . 

3.2 Gauge inveiriance at quadratic order 

In this subsection wc will derive the gauge invariant combinations to second order in the 
fluctuations. The idea is the same as in the previous subsection: we insert the harmonic 
expansion of the fluctuations and the gauge parameter into (3.7) (which now includes all 
terms) and read off the transformation of each coefficient. Then we seek a quadratic modi- 
fication of each field combination that is gauge invariant. One complication in this case is 
that because the r.h.s. of (3.7) is non-linear one needs to project onto the basic of spher- 
ical harmonics in order to extract the transformation of the coefficients. The analysis can 
be readily carried out in all generality but for the applications considered in this paper it 
is sufficient to consider only the modes that couple to scalar harmonics and their deriva- 
tives, i.e. we set to zero all fields that couple to vector and tensor harmonics (and their 
derivatives) (as well as C^^^)- Including these fields would result in additional terms in the 
gauge invariant combinations below. Since no confusion can arise we also drop the sub- 
script (s) from relevant fields and use the condensed notation for indices and arguments: 
^ (f)^,Ti'^'^ — Tr^,z{ki) — zi,a{ki,k2,ks) — 0123 etc. Note that we consistently use 
the notation ■0 to denote a field in the harmonic expansion of the supergravity field; ^ to 
denote a field which is gauge invariant to linear order and tp to denote the field which is 
gauge invariant to quadratic order. 
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3.2.1 Scalar fields 

We first discuss the scalar fields, ir^^ , (p^^^-^ and 6^^^. Their transformations are given by (we 
suppress the linear terms determined in the previous section and factors (C^^C^^C^^)): 

2'i>'edr,s + ih^e^' + e'D,^')a,2S + ie^' + 2e'Bl)h,s^ ; (3.16) 
(cV'ei23 + ie^Di,4>^ + ^eT^^)d2i3 + 2^^25^0123) ; (3.17) 

{ie'D^b^ + AH''e)ibu3 + A'ai23)) , (3.18) 

where the triple overlaps 0123 = a{ki, k2, ^a), &123 = ^'(^i) ^2, k^) etc are defined in appendix 
B.l and summation over {I2,H) is implicit. 

Prom these transformations one can infer quantities which are gauge invariant to quadratic 
order: 

1 ff2 2A^ \ f h} 

TT^ = TT^ - — M -A2ai23 + ^'123 " ^'^213) 4>'^T^^ + Ul23 " 2^^123+ (3-19) 



5tt^ = 








S(f>^ = 


1 




Sb^ = 


1 





A^( jA2ai23 + lbi23 - 
5 2 bqi 



,1 \ / ^ OA^ 

-d2i3) <f>''(f>^ + D^'7^^al23 + B^f^ib^s C123) 

Qi J ^ V H 

b^ = + ^ {^^'b%,2 + ^^d2n<P'-' + (^^'312 + 2^d2,3 + g^ei23) 0^0^ 

+Bl (^^'^ci23 + ^D%%,s^ ) . (3.20) 

The Ii = sector is special because the scalar harmonic is constant. Notice that we only 
have one scalar in this sector, namely tt^. Working out the gauge transformation yields 

JttO = z{k) (2^^^^qik) + ^A^^^Tf^ + e^D^TT^ - (e^TT^ + 2e^Bl)A'^ . (3.21) 

Prom here we obtain that the gauge invariant combination is (notice that tt'^ was gauge 
invariant to leading order, i.e. 7r° = 7f°) 

7r° = 7f° + zik) (^AV'vt^ - ^A^(A^ + 8)(l)'(t>' - B^^'D^t,' + A' Bi"' Bj)j , (3.22) 

where the summation over I is implied and tt^ and B>^^ are defined in (3.11) and (3.12). 

Let us now consider the field (^(j) associated with the tensor harmonic. Whilst this is 
gauge invariant to leading order, at the next order it transforms as 

41,) = ^ [e^P'e^L - ie'D,^' + le^' - 2e'^)c^^,^ , (3.23) 

where the normalization factor 2;(() and overlap integrals 0^23 ^tc are defined in the appendix 
B.l. Thus the gauge invariant combination to this order is 

<PU = Ht) + {{-B^'Bl + ivrV^' + Y^^'mc^, - ^ega^^^^) • (3-24) 
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3.2.2 Tensor fields 



We now turn to the non-scalar sector. As we will see in the next section the field equations 
algebraically relate the field bj^^p^ to the field 6^^^ (more precisely the field equations relate 
the corresponding gauge invariant combinations) so we need not discuss this field. Further- 
more, -B^j^^ is pure gauge. It is useful however to introduce the following combination that 
transforms nicely up to quadratic order 



= ^M + ^(-^(l^^/''^'^' + ^''^^-)«i23 (3.25) 



((l)^4>^{-hi2Z + ^ei23 + ^^213) + T7^d2130^7r3 + J( — Cl23 - ai2^)BlB''^\\ . 

V 8 8gi 5gi lOgi 2 qi J J 
This transforms as 

5Bl = il + ^ieBlbi23 + r'^.^3ai23). (3.26) 
Now consider the KK graviton modes, hjlj^. The gauge transformation reads 

6hl, = 1 [{D^e^hl + D^e^hl^ + e^Dji%)a^2Z + {.eh% + 2D^^eBl))hi2^) . (3.27) 

Prom this we derive the following gauge invariant combination (Ji ^ 0) 

K'^ = hl.-D^Bl-D,Bl-^^{^{<t>%, + ]^D^(t>''D,<t>'')bi2Z (3.28) 

+ {D^B>^^hl^ + D,B^%^ + B'^'^Djil, + D^B>^^D,Bl + B^^Hg'^, - BlB^aus) ■ 
Let us now discuss the /i = case. transforms as 

5hl, = D^e^'hl. + D^e^hl^ + e^Djil, (3.29) 

+ z{k) (D^e'hl + D,e'h{^ + e'Dx^u - + ^d^i^^'bi^)) ■ 

We introduce 

+D^B^'hi^ + D,B^%^ + B^'D^hl, + D^B^'D.Bi + B^'Big'^^, - B^B^) 

(the term linear in vr'' was added in anticipation of the fact that it is /i^j^ + \TT^g"^p that 
satisfies the linearized equations of motion, see the discussion around (4.5)). Recall that 
this mode is a correction to the spacetime metric 

9^tu = g^u + (3-31) 
so the combination should transform as 

= D^Q, + DIC^, (3.32) 



13 



where is the covariant derivative of the corrected metric (3.31). Indeed one finds this to 
hold with 

= + z{I){e'D^Bi - K'i'Bi). (3.33) 

4 Field equations 

The field equations for the coefficients in the harmonic expansion were derived in the de 
Donder gauge at linear order by [21] and at quadratic order in [26, 27, 28]. The gauge 
invariant variables derived in the previous section allow one to relax the gauge condition. 
Indeed notice that the gauge invariant variables when evaluated in the de Donder gauge 
become equal to the fields used in [21, 26, 27, 28]. It follows (and we have explicitly checked 
this in detail) that the field equations with no gauge condition imposed can be obtained 
from the results in [21, 26, 27, 28] by simply replacing each field with its gauge invariant 
generalization. 

4.1 Linear order 

In this subsection we summarize some of the results of [21] (a summary of the derivation 
is given in appendix C). As just mentioned, one can relax the gauge fixing condition by 
replacing all fields by the hatted versions given in the previous section. 
The scalars satisfy the following equations 

□s^i = k{k-A)s^\ k>2, 

□1^1 = (A; + 4) (A: + 8)^1, A;>0, (4.1) 

a^l^ = k{k + 4)4>ll^ k>2, 

where we introduce the combinations 

'" = 20(r+2)<*"-l"(' + ^)»")- «" = 20(r^2)(*"+^*»">- 

with inverse relations = -s^^ + f\ tt^^ = lOks^^ + 10{k + A)f^. 

The remaining modes that couple to scalar spherical harmonics are the KK gravitons. 
They are described by transverse and traceless fields 

(^fi - /ifi L. -D. D .(-TT-^i - 126^1), h 7^ 0. (4.3) 



satisfying the equation. 



iU-{k{k + A)-2))d].=Q, h^O. (4.4) 
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The 7i = case is special in that this mode describes a deformation of the background 
metric. The combination that satisfies the 5d hnearized Einstein equation is 

h',.' = {hi + (4.5) 

One can understand the origin of the shift by tt" by considering tlic reduction of the lOd 
action to five dimensions. Keeping terms hnear in the fluctuations, the volume of compact- 
ification manifold is 

J d^yVd^b = 7r3(l + IttO) (4.6) 
It follows that the reduced action is 

S5d - J Vd^t^((l + ^7rO)i? + • • •) (4.7) 

and a Weyl transformation is required to bring the action to the Einstein frame. The 
transformation from to h'^^^^^ is precisely this Weyl transformation. 

4.2 Quadratic order 

The derivation of the equations of motion to second order in fluctuations was discussed in 
[26, 27, 28] and is summarized in appendix C. For our applications it is sufficient to retain 
only the quadratic coupling to the field s^. 

4.2.1 Scalar fields 

The corrected field equation for the scalar fields ijj = {s^, s^, t'^, t^, t^, (/>^^^} is given by 

(□ - ml)ip' = D^22{^f + E^22D^s'Di'f + F^22D^^,D,^fD^i'D^)f, (4.8) 

where the coefficients D^22, E^22^F^22 can be obtained from the results in appendix C and 
are given in table 1. The fields entering the l.h.s of this equation arc the gauge invariant 
combinations to second order whilst the fields in the r.h.s. are the gauge invariant combi- 
nations to linear order (since the r.h.s. is quadratic in fluctuations). This follows from our 
general discussion and we have also explicitly checked that the terms involving (/)^^^ in the 
second order equations (with no gauge fixing imposed) are accounted for by the 0^^^ terms 
in the gauge invariant combinations. When is the leading non-zero field (as it is in the 
application discussed in this paper) the gauge invariant combinations take the form 

= V' + A^,^s^(t>l) + ^i>4'Ms)f + ^^sbD^'s'^bI)^ + A^bbB^^bI)^. (4.9) 
The coefficients A^sct>, A^4>4>, A^sB and A^bb are given in Table 1. 
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Table 1: Coefficients in (4.8), (4.9), (4.10) and (4.11). 



The field equations in (4.8) contain higher derivative terms on the r.h.s which can 
however be removed by the following transformation [26] 

* = w{^) + J^22{s^f + L^22Df,S^D^'s^) . (4.10) 

This transformation is the non-linear KK map to quadratic order in the fields. It maps 
solutions of the lOd fields equations to solutions of the 5d field equations, 

{n-ml)^ = X^22{S^f. (4.11) 

The coefficients w{il)), J ^22, L ^22 and Z^22 are given in Table 1. We include on the r.h.s. 
only the terms quadratic in S'^ because these are the terms that are relevant for us. We note 
however that all quadratic terms (and cubic scalar couplings [27]) have been determined in 
the literature [26, 27, 28]. The results in Tabic 1 arc in agreement with the results in these 
papers. The field equations can be integrated to a bd action, and the constants w{'ip) have 
been chosen such that the overall normalization agrees with the one in [8, 9] 

fj2 I- 1 1 

Sm = ^J d''xVG{-R + -G^^d^^d,^ + Vm. (4.12) 
Using the quadratic 10-dimensional supergravity action computed in [26, 29] one finds that 



, r, 8k(k-l)(k + 2)z(k) , h(t){k^*'^) , 
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= V (M) • 

From Table 1 we see that the only non-zero cubic term in the potential (that is also quadratic 
in S"^) is the cubic self coupling of S"^ and its coefficient, — 4/3-\/6, precisely agrees with 
the corresponding coupling in 5d gauged supergravity (after matching sign conventions), 
compare with (2.6) of [9]. 

It is important to note that the transformation (4.10) gives an explicit map between 
solutions of the ten dimensional equation and solutions of the five dimensional equation and 
vice versa, i.e. any solution of the five dimensional theory specified by the action (4.12) can 
be uplifted to a ten dimensional solution. We emphasize that this map is valid irrespectively 
of whether there is a consistent truncation since we keep all KK modes. 

4.2.2 Tensor fields 

Let us consider first the graviton. The quadratic correction to the gravitational equation is 
obtained in appendix C: 

{Le + ^)h% = ^ (^-'^D^DpD,s''D,DPD''s'' - ^{D^D^Dps'^DPs^ + D^DpfD^D'^s'') 

+^D^,fD,s^ + gl, [~BpB^s''DpB^%' + 2QDpS^DpS^ - ^s^^ ^ , (4.14) 

where Le is the linearized Einstein operator (C.20). This equation contains higher derivative 
interactions. Just as in the case of scalars one can remove them by considering the following 
transformation 

G,. = hl-^ (^iD.D^s'D^Dps' - ^{s^D.D^s' + (^(D^^)^ - ^{8^)9%) , (4-15) 

as can be verified using (C.21). In terms of these variables the field equation becomes 

R^,[G] = 2(T^, - ^G^.r;) (4.16) 

where is the Ricci tensor of G^i, and 

T^u = d^S'^d^S'' - G^4^idS''f + V{S'')). (4.17) 

The equation (4.16) is indeed the field equation for G^i, derived from (4.12) and T^^ is the 
corresponding matter stress energy tensor, where keeping with our approximations we only 
retain terms quadratic in S"^. 
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Now let us briefly describe the equations determining the KK gravitons. The traceless 
part of the ten-dimensional field is given by 

C) = + (^. + l)W3) °"-°-''^" - 

where all fields are the appropriate combinations which are gauge invariant to quadratic 
order. is transversal but i/jj^^^ is not; indeed its defining equation is 

L>'^D>^i = Z^H^], (4.19) 

where Z^^ [s] is quadratic in the field and is determined by the quadratic corrections to 
the trace of the Einstein equation in the AdS directions. The equation for the transversal 
field is then 

i^ + ^-Hk + mll.) = {'2LE + 8 + k{k + A))i;ll^^+Zll^y, (4.20) 



(2LE + 8 + A;(fc + 4))V^;f.\ 



where Zj^^^^^ is again quadratic in the field and follows from the corrections to the (fj^v) 
Einstein equation. is a transversal field which is quadratic in and contains up 

to six derivatives. We have verified that that the right hand side of the equation can be 
written in the latter form; this follows from the detailed structure of Z^^^^ and tpj^^. It is 
then immediately manifest that if one removes the higher derivative terms in the equation 
by defining the five dimensional field as ^jj^, = <;^^(^,^) — dimensional field 

satisfies the free field equation. This is in agreement with the result of [27] which found 
the corresponding cubic coupling to vanish and implies that the five-dimensional field must 
vanish to the order to which we work. As we discuss later, there is no physical content in 
these fields (to the order to which we work), so we suppress explicit details of the (rather 
complicated) KK reduction map. 



5 Holographic 1-point functions 
5 . 1 Generalit ies 

The KK reduction discussed in the previous sections provides an explicit map between ten 
dimensional solutions and five dimensional solutions as well as an associated five dimen- 
sional action for gravity coupled to massless and massive KK modes. If one would consider 
this problem in full generality the resulting action would involve an infinity of fields. For 
determining the holographic 1-point functions, however, we are only interested in the near 
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boundary behavior of the solutions (as is reviewed below). The near boundary expansion 

effectively decouples all but a finite number of fields, the number of which depends on the di- 
mension of the operator whose 1-point function one is computing; the higher the dimension, 
the greater the number of fields switched on. 

Starting from a five dimensional action there is a well developed method for comput- 
ing holographic 1-point functions, namely holographic renormalization [7, 9], see [17] for 
a review. Recall that the basic formula expressing the AdS/CFT correspondence [2, 3] 
relates the on-shcU supcrgravity action with prescribed boundary conditions for all fields to 
the generating functional of QFT correlators with the boundary fields playing the role of 
sources^, 

(exp (^-5qft[G'(o)] - J d^a;y^O(x)$(o)(x)^ = exp(-5sG[G(o), $(o)]) (5-1) 

where G(o) > ^(o) the fields parameterizing the boundary values of the bulk metric G and 
of other bulk fields denoted collecting by 

Naively, both sides of this relation diverge: the l.h.s. suffers from (the well known QFT) 
UV divergences and the r.h.s. suffers from IR divergences (due to the infinite volume of the 
spacetime). The divergences on the l.h.s may be dealt with by standard renormalization. 
The infinities on the r.h.s. are dealt by holographic renormalization. Namely, one adds a 
number of boundary counterterms that cancel all possible infinities that can arise in the 
on-shell action. Holographic 1-point functions in the presence of sources arc then obtained 
by computing in full generality the first variation of the renormalized on-shcU supcrgravity 
action. This leads to relations between the 1-point functions and certain coefficients in the 
near-boundary expansion of the bulk fields. This relation effectively replaces (5.1) since 
higher point functions can be computing by further differentiating the 1-point functions 
w.r.t. sources. 

The near-boundary expansion of the bulk metric G^i, and scalar field ^'^j where k is the 
dimension of the dual operator, take the form 

= ^ + ^ {G(o)ij{x) + z'^G^2)ijix) + z'^(G(^4)ij{x) + logz^/i(4)jj(x))) dx'dx^; 
^\x,z) = z2(logz%)(x) + 4.2o)(ar) + •••); 

$*^(x,z) = z(^-'=)$f„)(x) + ... + z%,_4)(x) + ---, k>2. (5.2) 

In these expressions the boundary fields G(o)y , 4*(o) ' ^(o) parametrize the Dirichlet bound- 
ary conditions and are also the field theory sources for the QFT stress energy tensor and 
■^We work with Euclidean signature. 
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operators of dimension 2 and k, respectively. The near-boundary analysis determines all 

coefficients in these expansions except the ones corresponding to the normalizable modes, 
namely G(4)ij,^fQ^:^^2k-4)- These are related to 1-point functions, as we review below. 

5.1.1 Radial Hamiltonian formalism 

The structure of the 1-point functions is most transparent in the radial Hamiltonian formal- 
ism [16, 10] (see [13, 14, 15] for earlier work). So before giving the explicit relation between 
the 1-point functions and the coefficients of the asymptotic solutions we digress to explain 
this relation. Let us define a radial canonical momentum for each field as 

dL , , 

where L is the Lagrangian and prime denotes differential w.r.t. r = — log z. A covariant 
version of the near boundary expansion in (5.2) is provided by the expansion of momenta 
in eigenfunctions of the dilatation operator, 

fe = £^A(2,«^+p-4)*'^) =-.^(l+0(.)), (5.4) 

where is the induced metric at the regulating surface z = t. The last equality follows 
from the leading asymptotics in (5.2). The near boundary expansion (5.2) now translates 
into the following expansions 

7r](x, e) = \/7(7r(o)j + • • • + vr(4)* + #(4)* log • • •); 

7r'=(x,e) = v^(7rf4_ifc) + --- + 4) + 4)loge2 + ...), k>2 (5.5) 

where vrj and vr^ are the radial momenta for the bulk metric and the scalar field , respec- 
tively, and each term in this expansion transforms as indicated by its index 

'^D7r(n) = -n7r(„), (5.6) 

except for 7r(4)*-, 7r^2)i ^(fe) which transform inhomogeneously, with the inhomogeneous term 
being equal to minus two times the coefficient of the logarithmic term, namely 

One advantage of the momenta expansion (5.5) over the near boundary expansion of the 
bulk fields (5.2) is that the momentum coefficients T^{n)ji''^\n) covariant w.r.t. 4d diffeo- 
morphisms that respect the regulating hypersurface z = e whereas the coefficients in (5.2) 
are not. 
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The coefficients in the momentum expansions can be obtained by inserting the expan- 
sions in ffamilton's equations. This leads to a number of equations obtained by collecting all 
terms with the same weight. One then solves these equations iteratively and each of them 
algebraically determines one of the coefficients in the expansion in terms of coefficients with 
lower weight. This determines all coefficients except vr(4)* and tt^^^-), just as the asymptotic 
analysis of the bulk equations determines all coefficients in (5.2) except for the normalizable 
modes. 

The renormalized 1-point functions are now simply given by the coefficient of the right 
dimension 

{Tij) = 7r(4)jj 

{O') = 4) (5.8) 

Following [16] one can show that there is a one to one correspondence between the momen- 
tum coefficients and the coefficients in (5.2). In particular, 

7^(2) = ^ (20(0) j , Trffc) = ^ ((2^^ - 4)0(2ik-4) + lower) (5.9) 

where the factor A/^^/27r^ is due to the overall factor in (4.12) and "lower" indicates terms 
with index less than {2k — 4). These terms are local functions of the sources so they are 
not important in computation of n-point functions for n > 1 (they lead to contact terms). 
They are important in the computation of vevs in cases where the solution describes a 
deformation flow [8, 20]. The specific example we discuss in this paper however is a vev 
flow so we need not specify them. 

5.2 bd supergravity fields 

The part of the bd action involving the metric and the field S"^ is same as the sector of 
gauged supergravity analyzed in [8, 9] (where S*^ was called <I>) (see also [15, 10]). Thus, 
the results for the 1-point functions carry over unchanged. The result for (O^) is as given 
above 

{O') = ^ (2^)) , (5.10) 

and for the stress energy tensor 

i'^ij) = 2^ {^{4:)ij + ^^{o)G{o)ij + g[Tr<^(2) - (TrG'(2))^]G(o)ij (5.11) 
~ l^i^l2))ij + 7G'(2)iiTrG(2) + -/i-(4)ii + ^^^fo) ~ ^fo)) ^(o)^ i_o)ij ■ 
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5.3 KK modes 

We now move to the one point functions of the other fields 5^, S^,T^, T'^,T^, Prom the 
last row of Table 1 we see that the cubic couplings (relevant for us) vanish for all of them 
so to the order to which we are working their field equations are just free field equations. 
One would therefore naively conclude that the one-point functions are simply given by 
(5.8)-(5.9). It turns out however that this is not true and there is an additional subtlety. 

The 1-point functions (5.8) were derived starting from a 5d action but in principle one 
should really evaluate on-shell the lOd action. In the majority of cases this distinction does 
not make any difference in practice, but there is a distinguished class of additional finite 
(non-local) boundary terms that one obtains from reducing the lOd action. These boundary 
terms are relevant for the computation of extremal correlators, namely n-point functions of 
1/2 BPS operators whose dimensions are {ki=k2 + ■ ■ ■ + kn, k2, ■ ■ ■ , kn} ■ These correlators 
are special in that they factorize into a product of 2-point functions. 

The bulk couplings associated with all extremal 3-point functions were shown to be zero 
in [26] (a result wc reproduced here for the coupling (S'^)^S"^, see Table 1). The 3-point 
functions of the corresponding operators however arc non-zero. It was shown in [22] in a 
specific example involving the dilaton and the t-field that even though the bulk contribution 
to the three point function vanishes, there are additional boundary contributions which 
lead to the correct 3-point function. By supersymmetry, the same should apply to all other 
extremal 3-point functions. It was further conjectured that these results generalize to all 
extremal n-point correlators. We refer to [30] for further discussion and references on this 
issue. 

These results should follow from holographic renormalization by starting from the lOd 
action, requiring that the variational problem is well posed and then KK reducing the 
action with boundary terms. Recall that in the 5d context all boundary terms, including 
counterterms, are uniquely fixed by the requirement that the variational problem is well 
posed with chosen boundary conditions [31]. We leave a detailed derivation for future work. 
Here we will instead use known results in order to fix the form of the 1-point functions. 

Covariance under 4d diffeomorphisms implies that the 1-point function should be a 
function of the coefficients in the momentum expansion (5.5). Furthermore the dimensions 
should match between the l.h.s. and the r.h.s. and the result for extremal ra-point functions 
should factorize into products of 2-point functions. This uniquely fixes the form of the 
1-point functions. For concreteness we discuss the 1-point functions of and but the 
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generalization is obvious. Thus 

{O') = 4)+a422(4))' (5.12) 
{O^) = vrfe) +a642vrfi)7r^2) +«633(7rf3))^ + 06222 (vr^2))^ (5-13) 

where 0422 , 0542 , oess > 05222 are numerical constants that we show how to compute in the 
next subsection. Let us explain the structure of these 1-point functions. The leading term 
TT^j,^ is the term discussed above in (5.8). Observe that the non- linear terms are possible 
only when the weight of the operator in l.h.s. can be written as a sum of weights of other 
operators, which is also the condition for extremal correlators. One could also consider 
adding terms involving tt^ with n < k. Such terms can only possibly contribute to the 
coincident limit of correlators (since tt^ with n < A; is local in the sources) or to vevs of 
solutions describing deformation flows, so they are not important for our analysis. 



5.4 Extremal couplings 



In this subsection we will compute the coefficients 04227 fi6222- The coefficients 0533 and 
0642 could be computed in similar way but we will not need their explicit values in this 
paper. To obtain the coefficients 0422, 00222 we compute the 3- and 4-point functions starting 
from (5.12) and then fix the coefficients so that the numerical factors agree with the the 
computation in free field theory. Note that the dependence on the coordinates is guaranteed 
to be correct by the structure of (5.12) (as should become clear shortly). 
We start with the computation of 0422- By definition 
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To evaluate these expressions we need to know tt^^^ to linear order in 0^^^ and vr^^^ to 
quadratic order in . Recall that tt^^^.^ is proportional to the vev part of the solution (see 
(5.9)), so in order to compute 2- and 3-point functions we need to solve bulk equations 
expanded around the background solution (which in the current context is just AdS) to 
linear and quadratic order, respectively, and then extract the coefficient of order z^. For 
extremal couplings the cubic coupling to is zero, A422 = 0, so the bulk equation for 
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continues to be □ 5^ = and the solution does not acquire dependence on 0^^^ so the second 
variation of vr^^^ w.r.t. 0^^^ is zcro^. Therefore only the last term in (5.15) contributes and 
using (5.14) we see that the extremal 3-point function is a product of 2-point functions 

{0\xi)0\x2)0''{x3)) = 2a422(C'(xi)02(x2))(02(xi)C»2(x3)) (5.16) 

Since this 3-point function does not renormalize one can compute it via free fields, which 
allows us to fix the proportionality constant 0422- 

In the large N limit (i.e. dropping non-planar contributions) the free field computations 
for the 2- and (extremal) 3-point functions yield 

0-W0-(,)) = ^^'—^^A—p, (5.17) 

/ ^ (27r)2'=l \xi - X2r \X1 - 3^3 1 ' 

where the operators are defined in (2.3) and in the extremal 3-point function ki = k2 + k^. 
The normalization factors A/fe are chosen such that the 2-point function of (5.17) agrees 
with the supergravity results; in particular given that for A; 7^ 2 the supergravity result is 



whilst for A; = 2 the result is 



<OWto)> = £(;^). (5.20, 



The normalizations are thus 



^2 = 2V2— . 

A 

Inserting (5.17) in (5.16) and comparing with (5.18) leads to 

(5.22) 

where in the last equality we have used the explicit value of the triple overlap (C4C2C2). 
This can be obtained from the formulae (B.13)-(B.14) by computing the overlap of y^^''^) 

and (y(2,o))2 

using the explicit expressions (B.9). 
^In the non-extremal case, the bulk equation reads = XkimS^S"^' and the r.h.s. induces a correction 
to S'' proportional to so the second variation of 77^*°' w.r.t. ^^q^ and ^^2) non-zero yielding 

the 3-point function, see [17] for a detailed discussion. 
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The computation of the coupling 0^222 is analogous. The bulk quartic coupling was 
shown to be zero in [32] and thus the only contribution to the extremal 4-point function 
comes from the last term in (5.13) 

4 

{0'^{xi)0\x2)0\x3)0\x4)) = 6ae4222 ll{0\xi)0\xk)) (5.23) 

k=2 

The free field result for extremal 4-point functions (in the planar limit) is 

■ (iTTj^'^l |Xi — X2I \Xl — Xs\ |Xl — X4I 

(5.24) 

where ki = k2 + + ^4, which fixes the proportionality constant to be 

where {C^C^C'^C'^) was computed using the following integral formula valid for extremal 
overlaps 



/ 



7^3 

(A;i + l)(A:i + 2)2*=!-!^ " ^ ' 

along with the explicit forms of the spherical harmonics. 

5.4.1 Summary 

To summarize we have shown that 1-point functions of the operators C>^,C^ are ^ 

//o6\ 6 I 4 2, / 3 \2 I -^6 3\/3 2 \3 

(C ) = 7r(e) + a642vr(4)7r(2) + a633(vr(3) j + j^-^{-k^^^) 

Notice that although we have used the 3- and 4-point functions on AdS to fix the couplings 
0422 and 06222, these 1-point functions hold for any solution of the bulk field equations. 

Notice also that these 1-point functions are compatible with the (conjectured) structure 
of near-extremal correlators. Recall that near-extremal correlators have weights k\ = k2 + 
■ ■ - kn — 2m with < m < n — 2. These correlators are conjectured (and checked through 
order g^) to be sums of terms each of which factors into products of lower-point correlators 
[33]. When m = 1 the correlator is called next-to-extremal and factorizes into a 3-point 
function and (n— 2) 2-point functions. One can easily check that this structure emerges 
from (5.27) after using the fact that the bulk coupling vanishes. (To be more precise, the 
®The non-linear terms in these relations may have the interpretation as operator mixing between single 
trace and multi-trace operators. 
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bulk coupling is known to vanish up to quartic order and is conjectured to vanish to all 
orders). For example the next-to-extremal correlator (O^O^O^O^) is given by 

{0\x,)0\x2)0\xsmx,)) = 2a422 f T1t4^%5V^) C^!kP\] +■■■ 

\'5'/'(0)(^2)¥(o)('T3) / y(Pfo){x4) J 

= 2a422 {{0\xi)0\x2)0\x3)){0\xi)0\x4)) + ■■■) (5.28) 
where the dots indicate permutation in X2,X3, x^. 

6 Coulomb branch solution 

6.1 Continuous distributions of D3 branes 

It is intuitively clear that the Coulomb branch of A/" = 4 SYM should be described by 
multi-center D3 brane solutions. Solutions describing N separated D3 branes solve the field 
equations that follow from the bulk supergravity action coupled to the worldvolume action 
of N (separated) D3-branes. It will be important for us to keep track of all normalizations 
factors so we set the stage by first reviewing some standard material. The bulk action is 
normalized as 

S=^j d'^xV^iR +•••), 2k' = {2n)\ar9l (6.1) 
and the worldvolume theory is given by 

^-source = ^J J d^a^ Cdbi{ct'')5{x'^ - X^(a'^)), (6.2) 

a=l 

where the Lagrangian for each D-brane is normalized as 

CDBiicr") = r3( Vdet(7 + 2WF) + •••), ^3 = 3^ , (6.3) 

where 7jj = d^X^ djX^ qmn is the induced metric and derivatives are with respect to the 
worldvolume coordinates a*. The D3-brane solutions take the form 

ds'^ = H{x±)-^/'dxf^+H{x±Y/'dxl (6.4) 
F = ^{dH-^ Auj\\-*±d±H) (6.5) 

where a;|| is the volume form in the worldvolume directions, *±_ and d± refer to the Hodge 
star and exterior derivative in the flat overall transverse directions and H is a harmonic 
function. We are interested in the case of a uniform distribution of N D3-branes on a two 
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dimensional disc of radius I. Approximating the distribution as a continuum distribution 

with density 

Pir) = ^0{l'-r% (6.6) 

where r is the radial coordinate in the plane of the distribution, the solution for the harmonic 
function is (see for example [34, 35]) 

where = A7TgsN{a')'^ , y arc coordinates in the four dimensional space transverse to the 
distribution of the D3 branes and f lies in the plane of the distribution. 
A change of coordinates 



y = rsm9, r = -s/P + r^ cos 9 (6-8) 



brings the solution into the form 

s2a / r4ji:2 



where 



ds'^ = -j^ f dxj + J + y {C'^de^ + sin^ Odn'^ + cos^ ^^^^ (6.9) 

F = L-"^ Ir^{l + ^sm^e)dr + -l'^r^sm2ed9\ Auii 

--1/ P _\ 

+ sin^ e cos 9— (X^(l + ^ sin^ 9)de - sin 29df A dJ^g A dcj) 
^4 y 2r^ 4r'^ / 



/2 _ 72 
= 1 + ^ sin^ 9, = 1 + ^ (6.10) 



Now note that if we rescale the four dimensional coordinates as 

(6.11) 

the metric has an overall factor whilst the five form has an overall factor. These 
factors combined with the prefactor of (6.1) result in the overall normalization factor of the 
five dimensional action (4.12); we can therefore suppress the L prefactors in the rest of this 
section. 

6.2 Asymptotic expansion 

We now wish to expand the metric near the boundary. A systematic way to do this is to use 
Gaussian normal coordinates centered at the boundary of AdS^ and then expand all fields 
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using the radial coordinate as a small parameter. This radial axial gauge can be reached 

by the charge of coordinates 

^ = z{l+a^z'' + a2Z^ + 0[zf) 

e = e + biz'^ + h2z'^ + 0[zf (6.12) 



where 



ai = ^{2 -We), a2 = ^cos^6'(5 + llcos26') 

1 5 
6i = — cos6'sin6', 62 = 29 (" sin 26" + sin 46*) (6.13) 



The metric then takes the form 

ds'' = ^ + L_^l + aiz'' + a2Z^)dxl (6.14) 
+d9^{l + Piz^ + P2Z^) + sin^ 9dnl{l + jiz^ + 72^^) + cos^ ^#^(1 + 6iz^ + (^2-2^) 

where the coefficients ai,a2 etc. depend on the angular coordinate 6. By scaling 

z^ zl (6.15) 

the leading metric becomes a unit radius AdS^ x 5^ and factors of I appear in the fluc- 
tuations. These factors can be easily reinstated in the final formulae so for simplicity we 
set Z = 1 for now. Using the explicit form of spherical harmonics in (B.9) and (B.18) the 
deviation of the metric from AdS^ x can be rewritten as in (3.4) with the following 

coefficients (valid up to terms of order z^), 

hUz) = ^^z%j, hl{z) = Vu(^-l + ^^z^y,„ h%{z) = -^-^z%„ 

^\z) = \z\ ^\z) = ^2{z--^lz% ^\z) = ^-^z\ 
Similarly, from the expansion of the five form we obtain 

4 _ 39V5 4 
\s) - —640"" 

There are several comments in order here. Firstly, the solution is not in the de Bonder 
gauge, as one can see from the fact that the scalar fields (\?^^^^ and (/>^^^ are non-zero. Sec- 
ondly, the expansion contains many more non-zero fields that one would naively expect. In 
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particular, there are non-zero KK gravitons, hfj and hfj (which are dual to the operators 

of the schematic form, TiF^F^X^ for k = 2,4 (sec Tabic 7 of [30])), scalar field (j)'^^-^ (that 
couples to a tensor harmonic and is dual to the operator TrXXXXX^) scalar fields 
(that are dual to TrF^F^X'' for k = 0,2,4) and scalar fields s^,s* (that arc dual to the 
operators TiX^ , k = 2,A). However, we know that in the CB flow only the operators TrX*^ 
get a vev. So what is the meaning of the values of the additional fields? 

To answer this question we should apply our map to obtain the gauge invariant five 
dimensional fields. As a first step we need to construct gauge invariant combinations. 
Using (3.19)-(3.20)-(3.22) and (3.24) and the definition of t'', we get 

t' = -r^^\ t' = -^z\ t' = -^z\ (6.18) 
128 ' 160 ' 160 ' ^ ^ 

2 \/3 2 \/3 4 4 37\/5 4 /2 A in\ 

s = — z z , s = z , (bf^-. = 0. (6.19) 

4 160 160 ^ ' 

The five-dimensional fields are obtained from these by the KK reduction formula (4.10) 
yielding 

ro = r2 = r^ = $2 0, = ^{z'' -^z^), 5^ = 0. (6.20) 

^ ' v6 6 

We thus see that that all additional scalar fields are equal to zero! The same is also true 

for the KK gravitons but we do not give the details here. We would like to emphasize, 

however, that a non-zero answer for these fields would not be a problem for the duality. 

The only cases where it would problematic is if the non-zero values correspond to a source 

or a vev. Note that all additional fields correspond to irrelevant operators and a non-zero 

source would not be consistent with AdS asymptotics. Furthermore, the corresponding vev 
part would appear at much higher power of z. To understand the (possibly non-zero) values 
note that the 5d equations are schematically of the form 

{a + ml)^ = X^22iS^? + ■ ■ ■ (6.21) 

where $ denotes collectively the fields other than , and the metric that are turned on. 
Any non-zero value for these fields would simply be induced by interaction terms - such 
non-zero fields just reflect the non-linear structure of gravity. In our case, it turns out that 
the couplings A$22 are zero, so the fields $ had to zero to this order. 

The fields that are important to understand at each order z'' are the ones which cor- 
respond to operators whose vevs can receive a contribution from the asymptotics at this 
order. In our case, these are the fields 5^, 5^ and the metric G^,^. The solution we discuss 
can be reduced to five dimensions using a "consistent truncation ansatz" (see [34, 36]). 
The reduced model involves the metric and S^. The expression for in (6.20) exactly 
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agrees with the asymptotic expansion of the 5c? solution, compare with (5.2) of [9] (and use 
$ = —S'^). We will return to the metric momentarily. 

The expression for S*^ in (6.20) is new information. The fact that it is zero comes out of 
non-trivial cancellations and at first sight is surprising since the vev of the dual operator is 
non-zero. From Tabic 1 we see that the coupling of to (S'^)^ is zero; this is an example 
of an extremal coupling 

nS^ = 0. (6.22) 

In this case however the vanishing of the coupling only explains the absence of logarithmic 
terms in the asymptotic expansion of S^. Logarithmic terms in the asymptotic expansion 
are related to conformal anomalies. Such conformal anomalies due to 3-point functions 
are possible when the couplings are extremal, see section 2 of [37]. They are however 
proportional to the sources so they evaluate to zero on the Coulomb branch, in agreement 
with the absence of logarithmic terms in the asymptotic expansion. In other words, the 
vanishing of the extremal couplings is required by the AdS/CFT correspondence and the 
structure of the conformal anomaly. Equation (6.22) allows for a homogeneous solution that 
is proportional to and one might have anticipated that the homogeneous term would be 
non-zero, since the vev of the dual operator is non-zero. We will resolve this issue in the 
next subsection. 

We now return to the spacetime metric. We see from (6.16) that the metric is corrected 
at the normalizable mode order. More precisely the combination that diagonalizes the field 
equations to linear order is h'^^ = (h^^ + ^TT^g'^^,) and in our explicit solution this is given 

by 

Kz = = (6-23) 

Naively this would imply that the dual stress energy tensor is non-zero. The solution, 
however, is supersymmetric so the vev of the stress energy tensor must be equal to zero. 

(We discussed this point earlier in the introduction.) As mentioned above the lOd solution 
can be reduced to five dimensions using a consistent truncation ansatz. The asymptotics 
of the 5d metric were given in the introduction in (1.3). As mentioned there, despite the 
non-zero coefficient of the term, the vev of Tij is zero because of additional contributions 
to the 1-point function. This does not immediately resolve the issue however because the 
metric in (6.23) does not agree with the metric in (1.3) (when both written in the same 
gauge)! 

The issue here is that /i^^ is not the correct 5d metric. Firstly, h'^^ does not transform 
correctly, i.e. as a five-dimensional metric. As derived in section 3.2.2 the combination 
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which transforms properly is h^^ in (3.30). Evaluating this formula for the case at hand 
gives 

hi = = ^z%,. (6.24) 

The five-dimensional metric is now obtained by using the non-linear KK map in (4.15). The 
resulting five-dimensional metric (including the background term) is given by 

This metric is not in the same gauge as the metric in (1.3). To correct for that we change 
coordinates as 2; = z{l + jxf^^^) ^^"^ metric becomes 

= ^ + Y2^^ - ^^')dxl (6.26) 

which precisely agrees with (1.3) (after reinstating the factors of 1)! Note that this coordinate 
change does not affect any other fields to the order we work to. 

6.3 Comparison with field theory 

Given the asymptotic expansions of the five dimensional fields we can now read off the 
vevs for the corresponding corresponding operators. These must agree with the field theory 
results discussed in section 2 because, as we explained there, A/" = 4 supersymmetry forbids 
any quantum corrections, so this computation is a test of the gravity/gauge theory duality. 

Using the formula (5.10) we see that the vev of the A = 2 operator is given in terms of 
the normalizable mode Sq (x) of the supergravity field as 

{<^')suGRA = ^2i^Si{x)) = -^/, (6.27) 

where we use the explicit result for our solution Sq = l/\/6 and we reinstated the factors 
of I. This 1-point function was previously derived in [8]. To compare with the field theory 
expression (set n = 1 in (2.11)), 

{0\,T = (6-28) 

we need to normalize operators and the action in the same way. The normalization of the 
operators was given in (5.21). Expanding the source action (6.2) to leading order in a' and 
taking into account the rescaling of the R^'^ coordinates in (6.11) the kinetic term for the 
scalar fields is normalized as 

T,L'J i'.(l(axn ={4^)jJ ^^(\m?). (6.29) 
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On the other hand, the field theory computation was done with canonically normahzed 

scalars. It follows that the radius a of the distribution in field theory (with canonically 
normalized scalars) is related to the radius I of distribution of D3 branes by 



a = yAz. (6.30) 

Using the normalization A/'2 in (5.21) and (6.30) we find precise agreement between the 
supergravity and field theory computations! 

Next consider the A = 4 operator. Given the result of the previous subsection that the 
normalizable mode of the corresponding supergravity field vanishes, the vev (5.27) receives 
contributions only from the term quadratic in 7r^2) • 

^ ^ V5A/|iV^^ 22^5 ^ ^ 

where we used (5.8) and (6.28). This is precisely the correct field theory vev! 



Now let us consider the A = 6 operator. The vanishing of Trf.x and irf^^ means that in 



this case the formula (5.27) reduces to just 



The latter of these terms evaluates to 

Given that in the A = 4 case, there was no contribution to the vev from the bulk super- 
gravity field one might have wondered whether the same was true in the A = 6 case, and 
indicative of a more general result. However, (6.33) only accounts for one fifth of the field 
theory result, and there must therefore be an additional contribution from the supergravity 
field dual to O^. To verify this one would have to extend our supergravity computations 
to one order higher, including quartic couplings, to extract the normalizable mode of the 
supergravity field S^. Note that the structure of the 1-point functions is such that the 
terms non-linear in momenta always give a contribution that is proportional to the QFT 
vev. It is a curious fact that up to at least (which is as far as we explicitly checked) the 
proportionality coefficient is a rational number, despite the fact that intermediate formulae 
contain square roots. 

The vev of the stress energy tensor was already computed in [8] (using (5.11)) and, 
as noted earlier, was found to be zero, in agreement with the fact that the solution is 
super symmetric. We have thus succeeded in showing that the vevs of all operators up to 
dimension A = 4 are correctly reproduced by supergravity! 



32 



7 Conclusions 



We have developed in this paper a systematic method for constructing the holographic map 
between the asymptotics of a ten dimensional solution and the 1-point functions of the dual 
QFT. Our main goal was to develop an unambiguous method that can, at least in principle, 
always be carried out. The main elements entering our construction are (i) the development 
of a gauge invariant version of KK reduction; (ii) construction of the KK map to non-linear 
order and (iii) application of holographic renormalization, including a proper treatment of 
extremal couplings. 

One-point functions can be derived rigorously starting from a 5d action via holographic 
renormalization. Our strategy for obtaining the 1-point functions dual to general KK fields 
was thus to reduce the field equations over the compact manifold and then use holographic 
renormalization. Recall that holographic renormalization relates the vcvs to coefficients in 
the asymptotic expansion of the 5d solution. So to compute the vevs starting from a lOd 
solution one has to understand quantitatively how the solution is reduced to five dimensions 
at the non-linear level. The point is that non-linear terms can give a contribution at exactly 
the same order (in a radial expansion) as linear terms. However, 1-point functions of 
operators of a given dimension can only receive contribution from non-linear terms involving 
fields dual to operators of lower dimension. 

The KK reduction map is constructed by first computing the fiuctuation equations 
around AdS^ x to a certain order in the fluctuation fields and then finding the field 
transformation that removes the higher derivative terms from these equations. This field 
transformation is the KK reduction map (to this order) and the resulting field equations 
are the 5d field equations. 

We would like to contrast our procedure with the procedure of "consistent truncuation" . 
In the latter one only keeps certain (typically low lying) modes in the KK reduction and 
then has to prove that the dynamics of these modes decouple from the rest. The existence 
of such a truncation is highly non-trivial and only holds for special compactifications. In the 
AdS/CFT correspondence consistent truncation maps to the closure of a subset of operators 
of the CFT under OPEs. In our discussion we keep all KK modes, so there is never an issue 
of consistency. We are however interested only in the asymptotic expansion of the resulting 
field equations. This effectively decouples all but a finite number of fields at each order in 
the expansion. In particular, for computing the vev of an operator of a given dimension 
only the fields dual to operators of the same or lower dimensions need to be kept. 

The need for a gauge invariant KK reduction stems from the fact that the KK reduction 
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is most efficiently done in a specific gauge, the de Bonder gauge, but in general explicit 
solutions will not be - and many known interesting solutions are not - in this gauge. Reaching 
this gauge would require finding a transformation that in general is not easy to obtain (at 
the non- linear level). Thus instead of gauge fixing the diffeomorphisms we construct gauge 
invariant variables. This allows us to immediately lift results derived in one gauge to 
another gauge: one can first obtain the KK map in the de Bonder gauge and then relax the 
gauge condition by simply replacing all fields by their gauge invariant generalization. The 
construction of the gauge invariant variables can be done systematically in the number of 
fields and we have done so up to second order in the fields. 

A final subtlety involves extremal couplings. Extremal correlators involve operators 
with the dimension of one of them equal to the sum of the dimensions of all the others. 
Such correlators are non-zero and are believed not to renormalize. A naive computation 
in supergravity however would give zero because the corresponding bulk coupling vanishes. 
It was argued in [22] that precisely in these cases there are additional boundary terms, 
originating from the higher derivative terms in the fluctuation equations, that one should 
take into account when evaluating the on-shell action and these yield the correct answer. In 
holographic renormalization one effectively replaces the on-shell action by renormalized 1- 
point functions in the presence of sources. These 1-point functions are valid for any solution 
of the field equations and higher n-point functions can be computed by further differentiating 
w.r.t. sources. Additional boundary terms in the 5d action, beyond the ones implied by 
the bulk 5d action via the variational problem, would manifest as additional contributions 
to the renormalized 1-point functions. The form of the additional terms is uniquely fixed 
by general principles. This leaves a few numerical coefficients to be determined and these 
can be easily computed by comparing the extremal correlators computed in weakly coupled 
M = 4: SYM and in supergravity. 

Combining these elements one obtains a well-defined holographic map. In our discussion 
we focused on solutions of IIB supergravity that involve only the metric and the self-dual 
five-form and asymptote to AdS^ x but the discussion readily generalizes to include all 
other fields or more generally to any theory with solutions that asymptote to AdSm x X, 
for some m and any compact manifold X. 

Let us summarize the steps involved in the construction of the map: 

1. Expand the solution (using a radial coordinate z as a small parameter) up to certain 
order and write the deviation from AdSm x X in terms of harmonics of X. 

2. If the solution is not in the de Bonder gauge, combine the fluctuations in gauge 
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invariant variables. 

3. Use the KK map to obtain the asymptotic expansion of the corresponding lower 
dimensional fields. 

4. Insert the coefficients of the asymptotic expansion in the renormalized 1-point func- 
tions to obtain the vevs. 

The asymptotic expansion of the lOd solution in general will contain non-zero terms for 
many coefficients. The only ones that carry physical information are the ones that have the 
correct leading radial behavior to correspond to normalizable or non-normalizable modes. 
The former give a contribution to the vev of the dual operator and the latter to the coef- 
ficient of the deformation of the QFT Lagrangian by the dual operator. We should note 
however that the remaining coefficients will generically contribute to vevs of higher dimen- 
sion operators via non-linear contributions to the holographic map. 

The higher the dimension of the dual operator, the higher the order needed in the 
gauge invariant variables and the KK map. So although this work solves the problem of 
computing the vevs in principle, in practice the method becomes cumbersome to carry 
out when sufficiently high dimension operators are involved (but since the procedure is 
algorithmic one could in principle computerize it). In this paper we explicitly worked out 
the map to the first non-trivial order. This is sufficient to compute vevs of operators up 
to dimension 4 and thus covers all relevant and marginal operators in four dimensions. As 
noted in the introduction, more efficient methods may be available when the solution has 
special properties. In this paper we mainly aimed at settling the issue of principle in full 
generality. 

To illustrate the general procedure we analyzed a solution that corresponds to a partic- 
ular point on the Coulomb branch of A/" = 4 SYM. This is an interesting example because 
the vevs are protected by supersymmetry and therefore the supergravity dual must repro- 
duce them exactly. The vevs corresponding to fields in gauged supergravity were previously 
computed in [8, 9]. Here we computed in addition the vev of the operator of dimension 4 
and exact agreement with the quantum field theory values was found! This constitutes the 
first non-trivial quantitative test of gravity/gauge theory duality away from the fixed point 
that involves a vev of an operator dual to a KK field. 

In the discussion so far we focused on how to compute vevs starting from a given 
lOd solution. The (inverse of the) holographic map can be used to see how spacctimc is 
reconstructed from QFT data. In particular, we see from our discussion that the vevs of 
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operators dual to KK modes provide a harmonic resolution of the compact space. Prom a 
more general viewpoint, notice that in the radial Hamiltonian formulation of holographic 
renormalization the vevs are associated with the radial canonical momenta conjugate to the 
sources. The holographic map therefore maps the field theory data to the phase space of 
the gravitational theory. It follows that these data are sufficient to uniquely determine the 
bulk solution, even though the explicit formulae only provide an asymptotic solution up to 
a certain order. 

One could thus "holographically engineer" duals of interesting quantum field theories 
by starting from the field theory vevs and using the holographic map. For this procedure 
to yield a smooth geometry, the vevs should clearly be large compared to the string scale. 
Even if this condition is satisfied, there is still no guarantee that a smooth geometry would 
emerge. For instance, it is well known that a necessary condition for a smooth supergravity 
dual is that the conformal anomaly of the theory at the UV fixed point should satisfy c = a 
in the large N and A limit [11]. Let us also note that even if a given theory has a smooth 
dual geometry, in practice it may not be easy to sum the asymptotic solution into this 
smooth solution. New tools that capture global issues of the correspondence may be needed 
to properly analyze this problem. It would be very interesting to explore this line of thought 
further. 

In many cases it is clear from the construction of the supergravity solution (with AdS 
asymptotics) what the corresponding gauge theory dual is. For instance this is the case if 
the solution is obtained via a near-horizon limit from another solution. Yet there are many 
other cases where solutions with AdS asymptotics have been obtained by directly solving 
the supergravity equations and there is no physical argument that would identify the dual 
theory. The work presented here should be useful both in verifying the gauge theory duals 
in cases where a proposed identification is available and also for extracting the gauge theory 
dual in other cases. 
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A Harmonic expansion of the antisymmetric tensor 



We expand the antisymmetric tensor as 

aabcd{x,y) = ^{bll^{x)eabc/DeY'{y) + bll^{x)eabc/Y,''{y)) (A.l) 
Gauge transformations act on the 4-form as follows 

SUMNPQ = '^D^mcnpq] 
The antisymmetric tensor parameter has the following expansion 

c^.p{x,y) = EcX(^)^'ny) (A.3) 



c 



fj.ua 



c,ab{x,y) = 

cabcix, y) = E(cg) {x)e^/^DaYt (y) + cj" (x)i?[„y/- {y)) (A.4) 
This implies the following gauge transformations for the fields, 

^Khp = Dpc[ly '^^w. = ^z'^" - K) = - Ml) 

It follows that the following combinations are gauge invariant, 

bjlupa = bXa + ^D^,Kl).pa] (A-6) 

(v)pup {v)pup 2 ''^ M^p] 

2 

uho _ uho I _ n, llio 
"{t)pv - "{t)pu ^ ^ l^^"{t)u] 

b{l)p = Kl)p+ {Ah-4)^i^Kh 

Indeed the field strength 

/mnpsr = ^D^MaNPRS] (A. 7) 

can be expressed in terms of these modes. 
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The gauge used in [21], 

D^'aaMNP = (A.8) 

amounts to setting to zero 

bj\ = b]\ = U]^ = h\\ = 0. (A.9) 

Our normalizations are such that the gauge invariant variables evaluated in this gauge agree 
with the parametrization in [21]. 

B Spherical harmonics 

The defining equations for the spherical harmonics are 

UyY^^ = A^iy^i, A^i = -A;(A: + 4), A; = 0,1,2,... (B.l) 

UyYt = h^'Yt, = _(fe2 + 4fc-l), fc=l,2,... 

Vlab) = ^'"^(i)' A'" = -(A:' + 4A:-2), fc = 2,3,... 

□,y[J,° ^ A''°Y^;°, A^io = -(A;2 + 4A;-2), A; = 1,2,... 



Useful identities for the scalar harmonics include 

A^ 

5 

ayD^aDb)Y' = {W + A')D^,D,)Y' 
OyDaY^ = (A^ + 4)L>„y^ 



D-D(aD,)Y' = 4(1 + — )D,y^ (B.2) 



B.l Spherical harmonics with 5*0(4) symmetry 

We introduce the following coordinates on 

ds^ = de^ + sin^ OdCll + cos^ edcp"^. (B.3) 

The differential equation (B.l) for the scalar harmonics is separable. Imposing 50(4) 

symmetry implies that the spherical harmonics depend only on 6 and (f). The general 
solution can then be expressed in terms of a hypergeometric functions, 

y(^'-)(e,<^) = C(„,^)y^(0)e^-'^ (B.4) 

where C(^n,m) is a normalization constant and the function y^(0) is given by 

yiix) = xH,F2(-^(fc - jmj), 2 + ^{k + jmj), 1 + jmj; x^) (B.5) 
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with X = cosd (there are also a second solution with leading behavior a;~l"*l but this 
solution does not reduces to a finite polynomial for any choice of the quantum numbers). 
The hypergeometric function reduces to a finite polynomial when either the first or second 
argument is zero or a negative integer. This leads to the following cases 

{k = 2l, m = 2n), {k = 2l + l, m = 2n + l) n e [-1,1], I e Z+ (B.6) 

with 

yin(^) = x2H^F2(-Z + |n|,2 + /+|n|,l + 2|n|;.T2) (B.7) 



y2i+\(^) = ^l^"+^lii^2(-/ + |n|,3 + / + |n|,2 + 2|n|; 



Particularly relevant for us here are harmonics that are also SO (2) symmetric which are 
given by 

-'^''"'•"•^ (s<-'" <"='"^") ■ 

The lowest harmonics are therefore 

Y^'''^ = ^(Scos^e-l), (B.9) 

y(4,o) ^ J_(i0cos^^-8cos2e + l), 
4V5 

y(6,o) = ^(35cos5^-45cos''0 + IScos^^- 1) 

8\/7 

The overall normalization in (B.8) has been chosen so that the harmonics are normalized 
as in [26], i.e. 

Recall that the scalar harmonics can be represented as 

y^i = C7f..,^x^i---x^'= (B.ll) 

where .x*" are Cartesian coordinates on and C^^...^^ is a totally symmetric traceless rank 
k tensor of 50(6). The normalization in (B.IO) corresponds to delta function normalization 
for the C^'s, i.e. 

(C^^C^^) = Cj...j^C^''i = S^^^\ (B.12) 
For the scalar harmonics we use the following definitions: 
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Jy^Y^W^-^ =a{h,k2,k3){C^'C^^C^--') (B.13) 

j D^^D^^Y^^DaY^^D^Y^^ = c{ki,k2,k3){C^'C^^C^^) 
j Y^^D^^D^W^^DaD^Y^^ = d{ki,k2,k3){C^'C^^C^^) 

J D^''D''^Y''{2DaDcY''D^^Df,)Y'' + DcY''DcD^^Dt,)Y'^) = e{k^MM){C''C'^C'^) 
where 

q{k) = A\4 + ^A^), a{k^, k2M) = 2) U ■ ^^'^"^^ 

5 (is + 2)122(^-2) Q;i!a2!a3! 

and T, = ki + k2 + ks, ai=^{k2 + k^ — ki) etc. One can derive explicit formulae that express 
b{ki, k2, fcs), c(A;i, ^2,^:3), d{ki, k2, ks), e{ki, k2, ks) in terms of a{ki, k2, ks) by use partial in- 
tegrations. Useful identities include: 

d{k2,ki,k3) + c{ki,k2,k3) + ^^^b{k2,ki,k3) = 0; 
b{k2,ki,k3) + b{ki,k2,k3) + A^^a{ki,k2,k3) = 0. 



We further have 

D^qDq^Y^ = ^(1 -2 cos 29), D^qDq^Y^ = ^{2 + cos2e -5cos4e). (B.15) 
For the tensor harmonics we use the following definitions 

J Y^^DaY'^D.Y'^ = c^'\k?,k2, k3); (B.16) 
J Y^l{2DaD,Y''W(,Dh)Y'^ + D,Y'^D^D^^D,,)Y'-') = e^'\kf\k2,k3). 
The normalization of the spherical harmonic is defined as 

The only tensor harmonic of relevance here is 

yS'""^ = -3, Y^l'"^ = cos^ e{-3 + 15 cos^ 9), Y^^l = sin^ 9{2 - 5 cos^ 9) (B.18) 

where tp"- are the coordinates on S^; we thus do not need to discuss the tensor harmonics 
more generally. 
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C Field equations up to second order 



We discuss in this appendix the derivation of the field equations up to second order in 
fluctuations. 

The linearized equations read 

_ I '^iKL-po rpo M1M2M3 fn -\^ 

^MN = ^MN "I" ^MKMiM2M3-'^NL {^-^J 

^ tf jpO M1M2M3M4 . f jpO MlM2M3M4\ r. 

— -^[jMMiM2M3M4-t'N + jNMiM2M3M4^ M j=U 

^\Ji...Mz = ~ f*)Mi...M5 + 2^i-^Mi...M6 - ^h^Mi^M2...M5]K = ^ (^-2) 

where 

^MN — ^KhMN ~ -^DMDNh^-, ^MAT = 2^^^^^ ^Nh^ — D^Hmn)- (C3) 

Projecting these equations onto the various harmonics leads to the following equations^ 

{eS)\d,^d,,y'.=0 => {\k'' + ^ic^^) = 0; (C.4) 

(^i?)ly/i4=0 =^ i(D + A^"-2)<^5;=0; 

-LU O 

These equations lead to the scalar field equations quoted in section 4. 1 upon elimination of 
bixupa and h'^ and then diagonalization. 

We now move to the quadratic order. The field equations are 

_ t(2) _ 7.(2) 

where the quadratic corrections are given by [27] 

^m!...Mb = -2^LfMi...M5 + ^h[Mj*M2...M5\K i^-^) 
^ L K ^ L 2 ^ ML K K 

T^(2) _ _ r{2) 7\ 

, '^LXLrSpo po M1M2M3 ."^^rKiSi 1^X2821^0 po M1M2 

+ l^L^ MNMiM2M3^ NS "3]""' ^ MKiK2MiM2-'^ NS1S2 

^i,KS(f jpo M1M2M3 I J? TTio MiM2M3\ , 1 /• J? M1...M4 

— gj'^ {JMKMiM2M3-fNS + jNKMiM2M3^ MS ) + ■^JMMi...M4jN 



"^In comparing with [21] one should note that we expand in harmonics hfn, rather than h'^^ (compare 
(2.5)-(2.7) with our (3.4)). 
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where 

Rmn = -Dxih^hj^r^) + ^i^iv(/ii^L^M/i^^) + ^/^^iv^i^^i - ^W^^l (CS) 

These quantities were computed to second order in the field s in the de Bonder gauge in 
[26] , by substituting the hnear solution of the field equations 

Ki = U{k)s'^g% + W{k)D^^D,)s'^- (C.9) 
V{k) = -lu{k) = 2k, W{k) = ^ 



3 ' ' ' ^ ' fe + l 

As discussed in the main text, the resulting field equations will be applicable to other gauge 
choices provided that one replaces each field by the corresponding gauge invariant field. 
In particular, for the field "s^^" must denote the appropriate gauge invariant field. For 
computing the quadratic corrections, however, it is sufficient to use the field s^^ which is 
gauge invariant to linear order, since the difference between this field and the gauge invariant 
field is itself quadratic in fluctuations. 

C.l Scaleir fields 

To compute the corrected scalar equations we will need to use the following components of 
(C.6) and (C.7): 

Q\ ^ ^(^;?)lD(.D,)yi; (c.io) 

= 7^77^ ((C123 + ^231 + d32l)T23 + ^20x2^^ S^) , 

T23 = i3V2V3 + 5U2U3)sH^ + W2W3D^f'D''h'^D^^D^)s\ 

The notation in the first line implies the projection of the tensor (which is quadratic in 

spherical harmonics) onto the spherical harmonic. (Note that the factor of five is included 
in the definition so as to match the conventions of [26]). Similarly one has 

Ql = ^{tP «)|yi = ^ (I05i23 + T23(6i23 - 2/30123) + S2Di,s^ D^' s^us) , 
S123 = VsU2ai23D^'{s^D^s^) + W2Vsai23D^{D^f'D''h^D,s^) (C.ll) 

-V2V3b213sH^ - 8(ai23/2/3s's3 + ^230^^ D^s"^) - 0123^^2(64/3 + SOFs)^'^^. 

We also need 
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Ql^ = (iU2 + 3V2)s'^D^s^ + W2D^^Dfh'^DpsA 6213; 

(^/l^p(TT)|yi = Qhl^i^par, (C.13) 
Ql = -^(r23-(16F2/3 + 40y2V^3)s's^)ai23- 

For the scalar field coupling to the tensor harmonic we need 

Q'(l) = (7^i?)ly^i4, (C.14) 

^ ' ab 



Then the corrected equations of motion are written in terms of quantities just defined 

as 

(□_fc(fc-4))/i = -^l—^((k + A)ik + 5)Qi + Q2 + {k + A)iD^Q^ + Q,)y'- 
{a-{k + 8)kik + A)y^ = -^l-^{kik-l)Qi + Q2-k{D^Q^ + Q^)Y\ (C.15) 

whilst the corrected equation for the scalar (j)^^^^ is 
C.2 Tensor fields 



{a-k{k + 4M},'=2Q'}^. (C.16) 



(2) 

To compute the correction to the metric and KK gravitons we we need Tjli, . For brevity 

we will include only the terms of interest here, namely s^. The curvature contribution to 

iT^V)\Yi is then 

1/4 32 8 

- {--ai22D^DpD^s'^D^DPD''f - —ai22DPfDpD^D^f - -hi22D^Dpf D^DPf 

+^ai22{DPD^s^DpD^s^)g"^, + (^6/22 - S2ai22)s^D^D,s'' (C.17) 



whilst the field strength contribution to )|y/ is 

^ (^^.(f a/22(2(s2)' - DpD,s''DPD''f) - 4bi22{DpS^DPs^)) - ^ai22§''D^D,s'' 
+8D^s'^D^s%22) ■ (C.18) 

These lead to the following equation for the graviton, 

{Le + A)h% = T^J \yo - g%{\Ql + 8QI)\yo (C.19) 
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where the hnearized Einstein operator is defined as usual by 

LeK. = ^(-DA^. + DpD^XP^ + DpD.XP^ - D^DA'^). (C.20) 

The term proportional to Q2 and Q4 arise when eliminating Dtt^ and h'^^ from the equation. 
The following identities prove useful: 

LE{{DPD''fDpD^f)gl,) = ~U{DPD'' s^DpD^f)g';,, - iD^D^DpD^fOP D'' s"" 

LE{D^DpS^D^DPs^) = 2DPD^s'^DpD^s'^ - 0^,0^0" fO^DpD^s'^ 

+g%{DPD''s^DpD^s^) - QDPs^D^D^DpS^ - ID^s^D^f 
+m^DpD,s^ - {DPs^DpS^)g%, (C.21) 
LE{{DPs^Dps'')g;,) = S{DPfDps'')g%- g^iDPD^s^DpD^f) 

-?,DPfD^D^Dpf - ?,Di'D^fDpD^f- 
Le{^Dp,D,s^) = -SD^s'^D.s'' - {Dfs^Dps'')gl, + D^D^s^D^DpS^; 
LE{{^?g%) = ^{^?g% - {Dn''Dpf)g% - ^s'^D^D.s'' - ZD^fD,s\ 
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